Fast Non-Adiabatic Two Qubit Gates for the Kane Quantum Computer 
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In this paper we apply the canonical decomposition of two qubit unitaries to find pulse schemes 
to control the proposed Kane quantum computer. We explicitly find pulse sequences for the CNOT, 
swap, square root of swap and controlled Z rotations. We analyze the speed and fidelity of these 
gates, both of which compare favorably to existing schemes. The pulse sequences presented in this 
paper are theoretically faster, higher fidelity, and simpler. Any two qubit gate may be easily found 
and implemented using similar pulse sequences. Numerical simulation is used to verify the accuracy 
of each pulse scheme. 



I. INTRODUCTION 



The advent of quantum algorithms [jj, |2j that can out- 
perform the best known classical algorithms has inspired 
many different proposals for a practical quantum com- 
puter @, 0, S IE II 13- One of the most promising 
proposals was presented by Kane 0. In this proposal a 
solid state quantum computer based on the nuclear spins 
of 31 P atoms was suggested. Although initially difficult 
to fabricate, this scheme has several advantages over rival 
schemes 0:0, ID, 00 Hi- These include the comparatively 
long decoherence times of the 31 P nuclear and electron 
spins 

[HIIllHEQHEinithe similarity to ex- 
isting Si fabrication technology, and the ability to scale. 

There have been two main proposals for pulse se- 
quences to implement a CNOT gate on the Kane quan- 
tum computer. In the initial proposal an adiabatic 
CNOT gate was suggested. Since that time the de- 
tails of this gat e have been investigated and optimized 
IB fisL I20L l2lj . This adiabatic scheme takes a total 
time of approximately 26/xs and has a systematic error 
of approximately 5 x 10 -5 [l^. As good as these results 
arc, non-adiabatic gates have the potential to be faster 
with higher fidelity and allow advanced techniques such 
as composite rotations and modified RF pulses [22I l23l | . 

Wellard et al. |24| proposed a non-adiabatic pulse 
scheme for the CNOT and swap gates. They present 
a CNOT gate that takes a total time of approximately 
80/xs with an error (as defined later in Eq. (jHU) of ap- 
proximately 4x 10 -4 . Although this gate is non-adiabatic 
it is slower than its adiabatic counterpart. For the non- 
adiabatic swap gate a total time was calculated of 192/is. 

One of the most useful tools in considering two qubit 
unitary interactions is the canonical decomposition 
[2^, [2^. This decomposition expresses any two qubit gate 
as a product of single qubit rotations and a simple inter- 
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action content. The interaction content can be expressed 
using just three parameters. In the limit that single qubit 
rotations take negligible time (in comparison to the speed 
of interaction), this decomposition can be used to find 
optimal schemes |2(1 l27j . and of particular inspiration 
to this paper is an almost optimal systematic method to 
construct the CNOT gate [28| . 

It is not possible to apply those optimal schemes 
|2fl l27j | directly to the Kane quantum computing archi- 
tecture. They assume single qubit gates take negligible 
time in comparison with two qubit interactions, whereas 
on the Kane architecture, they do not. Secondly, in the 
proposal for the Kane computer, adjacent nuclei are cou- 
pled via the exchange and hyperfine interactions through 
the electrons, rather than directly, and so we have a four 
'qubit' system (two electrons and two nuclei) rather than 
a two qubit system. Although we cannot apply optimal 
schemes directly, in this paper we use the canonical de- 
composition to simplify two qubit gate design. 

Apart from being simple to design and understand, 
gates described in this paper have many desirable fea- 
tures. Some features of these gates include: 

1. They are simpler, higher fidelity and faster than 
existing proposals. 

2. They do not require sophisticated pulse shapes, 
such as are envisioned in the adiabatic scheme, to 
implement. 

3. Any two qubit gate can be implemented directly 
using similar schemes. This allows us to implement 
gates directly rather than as a series of CNOT gates 
and single qubit rotations. 

This paper is organized as follows. Sec. [H] gives an 
overview of the Kane quantum computer architecture 
and single qubit rotations. Sec. IIIII describes the canon- 
ical decomposition as it applies to the Kane quantum 
computer. Section llVl describes pulse schemes for Con- 
trol Z gates and CNOT gates. Sec. |V| gives potential 
pulse schemes for swap and square root of swap gates. 
Finally, the conclusion, Section IVT1 summarizes the find- 
ings of this paper. 
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II. THE KANE QUANTUM COMPUTER 
A. The Kane Architecture 

A schematic diagram of the Kane quantum computer 
architecture is shown in Fig. ^ The short description 
given here follows Goan and Milburn ^8|. This archi- 
tecture consists of 31 P atoms doped in a purified 28 Si 
(/ = 0) host. Each P atom has nuclear spin of / = |. 
Electrodes placed directly above each P atom are referred 
to as A-Gates, and those between atoms are referred to as 
J-Gates. An oxide barrier separates the electrodes from 
the P doped Si. 



Under typical operating conditions, a constant mag- 
netic field B will be applied to the entire system, perpen- 
dicular to the surface. This contributes Zeeman energies 
to the Hamiltonian: 

Hb = —g n /J. n BZ n + hbBZ p . (3) 

A typical value for the Kane quantum computer of B = 
2.0T gives Zeeman energy for the electrons of hbB « 
0.116meV, and for the nucleus g n [i n B ps 7.1 x lCP 5 meV. 

The hyperfine interaction couples between nuclear and 
electronic spin. The contribution of the hyperfine inter- 
action to the Hamiltonian is 
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Ha = Aa e ■ a ni 



(4) 



where strength, A, of the hyperfine interaction is propor- 
tional to the value of the electron wave-function evalu- 
ated at the nucleus, 
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FIG. 1: The Kane Quantum Computer Architecture 



Each P atom has five valence electrons. As a first 
approximation, four of these electrons form covalent 
bonds to neighboring Si atoms, with the fifth forming 
a hydrogen-like S-orbital around each P + ion. This elec- 
tron is loosely bound to the P donor and has a Bohr 
radius of a* B ~ 3nm, allowing an electron mediated in- 
teraction between neighboring nuclei. 

In this paper nuclear spin states will be represented 
by the states |1) and |0). Electronic spin states will be 
represented by | |) an d | !)■ Where electronic states are 
omitted, it is assumed that they are polarized in the | [) 
state. X, Y and Z are the Pauli matrices operating on 
electron and nuclear spins. That is 



X 



,Y = <Jy,Z 



ID 



Operations which may be performed on any system 
are governed by the Hamiltonian of the system. We now 
describe the effective spin Hamiltonian for two adjacent 
qubits of the Kane quantum computer and give a short 
physical motivation for each term which makes up the 
overall Hamiltonian: 



A typical strength for the hyperfine interaction is A = 
1.2 x 10 -4 meV. Charged A-Gates placed directly above 
each P nucleus distort the shape of the electronic wave- 
function thereby reducing the strength of the hyperfine 
coupling. The nature of this effect is under numerical in- 
vestigation |29j . For the purposes of this paper we have 
assumed that it will be possible to vary the hyperfine 
coupling by up to approximately 50%. 

The exchange interaction couples adjacent electrons. 
Its contribution to the Hamiltonian is: 



Hi 



J<?e. 



(6) 



where e\ and e2 are two adjacent electrons. The magni- 
tude, J, of the exchange interaction depends on the over- 
lap of adjacent electronic wave-functions. J-Gates placed 
between nuclei distort both electronic wavefunctions to 
increase or decrease the magnitude of this interaction. A 
typical value for the exchange energy is 4 J = 0.124meV, 
and for the purposes of this paper we assume that it will 
be possible to vary the magnitude of the exchange inter- 
action from J = to J w 0.043mcV. 

A rotating magnetic field, of strength B ac rotating at 
a frequency of tu ac can be applied, perpendicular to the 
constant magnetic field, B. The contribution of the ro- 
tating magnetic field to the Hamiltonian is: 



2 

H = J2h Bz + H Af + Hj + H aCi . (2) 



H ac = -gn^nBac [X n COs(uJ ac t) + Y n SUl(w ac t)] 

+Mb 5 qc [X e cos(w ac t) + Y e sin(a; ac t)] , (7) 



where the summation is over each donor atom in the where the strength of the rotating magnetic field is envi- 
system, i. sioned to be B ac w 0.0025T. 
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At an operating temperature of T — lOOmK, the elec- 
trons are almost all polarized by the magnetic field. That 
is: 



-j « 2.14 x 10" 



(8) 



We assume that electrons are polarized in the | J.) state, 
and use nuclear spin states as our computational basis. 



B. Z Rotations 

Single qubit rotations are required to implement the 
two qubit gates described in this paper, as well as being 
essential for universality. In fact, as we will see they 
contribute significantly to the overall time and fidelity of 
each two qubit gate. It is therefore important to consider 
the time required to implement Z , X and Y rotations. 

In this subsection we describe how fast Z rotations may 
be performed varying the voltage on the A-Gates only. 
A Z rotations is described by the equation: 



(9) 



A Z gate (phase flip), may be implemented as a rotation. 
It is given up to a global phase by: 



Z 



-iRz (tt) 



(10) 



Under the influence of a constant magnetic field, B, to 
second order in A [Isj . each nuclei will undergo Larmor 
precession around the Z axis, at frequency of 



hwi = 2g n fi n B + 2A 



2A 2 



HbB + g n /j, n B 



:id 



Z rotations may be performed by variation of the hy- 
pcrfine interaction from A to A z giving a difference in 
rotation frequency of 



2{A-A Z ) + 



2(A 2 - Ap 
HbB + g n ii n B' 



(12) 



Perturbing the hyperfine interaction for one of the 
atoms, and allowing free evolution will rotate this atom 
with respect to the rotation of the unperturbed atoms. 
The speed of single atom Z rotations depends how much 
it is possible to vary the strength of the hyperfine inter- 
action, A. For numerical simulation we use the typical 
values shown in Table [I] 

Under these conditions a Z gate may be performed on 
a single nuclear spin in approximately 



t z « 0.021/is. 



Description 


Term 


Value 


Unperturbed Hyperfine 
Interaction 


A 


0.1211 x 10" 3 meV 


Hyperfine Interaction 
During Z Rotation 


A z 


0.0606 x 10~ 3 meV 



TABLE I: Typical Parameters for a Z Rotation 



These rotations occur in a rotating frame, that pro- 
cesses around the Z axis with a frequency equal to the 
Larmor frequency. We may have to allow a small time 
of free evolution until nuclei that are not affected by the 
Z rotation orientate themselves to their original phase. 
The time required for this operation is less than 



t F < 0.02^s 



X and Y Rotations 



(14) 



In this section we show how techniques, similar to those 
used in NMR 0, |3(], OOI > mav be used to implement X 
and Y rotations. X and Y rotations are described by the 
equations: 



Rx(e) 

Ry{9) 



e" 2 



x 



(15) 
(16) 



X and Y rotations are performed by application of a 
rotating magnetic field, B ac . The rotating magnetic field 
is resonant with the Larmor precession frequency given 
in Eq. (fTTfl . that is: 



(17) 



In contrast to NMR, in the Kane proposal we have 
direct control over the Larmor frequency of each individ- 
ual P nucleus. By reducing the hyperfine coupling for the 
atom we wish to target from A to A x we may apply an 
oscillating magnetic field that is only resonant with the 
Larmor frequency of only one of the atoms. This allows 
us to induce an X or Y rotation on an individual atom. 
To the first order, the frequency of this rotation may be 
approximated by: 



'iWi = 9n^nB ac 1 



gnl^nB 



(18) 



The speed of an X rotation is directly proportional to 
the strength of the rotating magnetic field, B ac - As the 
strength of the rotating magnetic field, B ac increases, the 
fidelity of the operation decreases. The reason is that in 
frequency space the Full Width Half Maximum (FWHM) 
of the transition excited by the rotating magnetic field in- 
(13) creases in proportional to B ac . That is, as B ac increases 
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Description 


Term 


Value 


Unperturbed Hyperfine 
Interaction 


A 


0.1211 x 10~ 3 meV 


Hyperfine Interaction 
during X Rotation 


A x 


0.0606 x 10~ 3 meV 


Constant Magnetic Field 
Strength 


B 


2.000T 


Rotating Magnetic Field 
Strength 


B ac 


0.0025T 



TABLE II: Typical Parameters for an X Rotation 



we begin to excite non-resonant transitions. The larger 
separation, in frequency space, between Larmor frequen- 
cies, the smaller this systematic error. Since the Larmor 
precession frequency depends on how much we are able 
to vary the hyperfine interaction, A, it determines how 
strong we arc able to make B ac 

For the purpose of simulation, the typical values shown 
in Table [H] for the unperturbed hyperfine interaction 
strength A, the hyperfine interaction strength during the 
X rotation A x , applied magnetic field strength B, and 
rotating magnetic field strength B ac were used. 

Using these parameters this gives the overall time to 
perform an X gate on a single qubit in approximately: 



tx ~ 



(19) 



Any single qubit gate may be expressed as a product 
of X, Y and Z rotations. Ideally, X and Y rotations 
should be minimized because Z rotations may be per- 
formed much faster than X or Y rotations. For example, 
a Hadamard gate may be expressed as a product of Z 
and X rotations: 



*-M5)MI)MI 



(20) 



Thus, from the above discussion, Hadamard gate takes a 
time of approximately: 



tn ~ 3.2/xs. 



D. Nuclear Spin Interaction 



(21) 



In this section wc show the results of second order per- 
turbation theory to describe the interaction between two 
neighboring P atoms. This interaction between nuclei is 
coupled by electron interactions. We consider the case 
where the hyperfine couplings, between each nucleus and 
its electron are equal, that is: 

A = A 1 =A 2 . (22) 

We allow coupling between electrons, that is: 

J > 0, (23) 



but restrict ourselves to be far from an electronic energy 
level crossing, 



(24) 



Under these conditions electrons will remain in the po- 
larized | II) ground state. 

In this situation analysis has been performed using sec- 
ond order perturbation theory [l^. To second order in 
A, the energy levels are: 



Eli 



in) 



E 



K) 



E, 



loo) 



-2fi B B + J + 2g n fx n B + 2A, 
OA 2 

-2fi B B + J — -, 

\x B B + g n \x n B 

2A 2 

-2ix B B + J — — — , 

fi B B + g n H n B - 2 J 

-2fi B B + J - 2g n fi n B - 2A 

2A 2 2A 2 



li B B + g n [i n B - 2.1 [i B B + g n HnB 



(25) 
(26) 

(27) 
(28) 



where the symmetric \s n ) and anti-symmetric \a n ) energy 
cigenstates are given by: 



K) = -J=(|io> + |oi», 
K) = i(|io>-|oi». 

Notice that the energies arc symmetric around 

A 2 

E a = -2fi B B +J 



(29) 
(30) 



\x B B + g n [i n B - 2 J 
A 2 



fx B B + g n fx n B 



(31) 



Since we are free to choose our zero point energy to be Eq 
(or equivalently ignore a global phase of a wavefunction, 
|^>)) we may rewrite the second order approximation as: 



E 

E\ 
E 



llDlii) 
I1DK) 

II) K) 

III) 100) 



Huj b , 
huj s , 
huj s , 

huj Bl 



(32) 
(33) 
(34) 
(35) 



where lu b and los are given by: 

hijj B = 2A + 2g n /i n B + 
A 2 



A 2 



fx B B + g n fJ, n B n B B + g n fi n B - 2,1 

A 2 A 2 

HbB + g n p n B - 2 J f-L B B + g n [i n B 



,(36) 
.(37) 



The reason for this representation of the energy will be- 
come clear in the next section. Typical values were used 
during numerical simulation of the interaction between 
nuclei are shown in Tabic ITTT1 
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Description 


Term 


Value 


Hyperfine Interaction 
during Interaction 


A v 


0.1197 x 10" 3 meV 


Exchange Interaction 
during Interaction 


Ju 


0.0423 meV 



then a similar analysis shows that 



TABLE III: Typical Parameters during Interaction 

III. THE CANONICAL DECOMPOSITION 

In this section we describe the canonical decomposi- 
tion, and describe how this decomposition may be ap- 
plied to the Kane quantum computer. 



Mathematical Description of Canonical 
Decomposition 



The canonical decomposition |25l |27| decomposes any 
two qubit unitary operator into a product of four single 
qubit unitaries and one entangling unitary. 



U={V 1 ®V 2 ) U can {Wx®W 2 ). 



(38) 



where Vi, V2, Wi and W2 are single qubit unitaries, and 
U can is the two qubit interaction. The symbol ® repre- 
sents the tensor product of two matrices. 

U can has a simple form involving only three parame- 



ters, a x , a y and a z 



JciyYtgiY ia z Z<SiZ 



(39) 



This purely non-local term is known as the interaction 
content of the gate. It is not difficult to show that each of 
the terms in the interaction content, e la * x ® x ; ^a y Y®Y 
and e !tt2Z8Z , commute with each other. 

Physically each of the terms e la * x ® x , e la y Y ® Y , and 
e ia z z®z correS p 01K i to a typ e of controlled rotation. For 
example, following Ref. 



e 



ia z Z<g>Z 



= cos a z I (£) I + i sma z Z £g) Z 
= cosa z (|0)(0| + |l)(l|)®7 

+isina 2 (|0)(0| - |1)(1|)®Z 
= |0)(0|®e M ^ + |l)(l|®e- M ^ 
= (I®e ia * z ) (|0)(0|®7 + |l)(l|®e- j2Q ^)(40) 

This shows that up to a single qubit rotation, e m,28Z 
is equivalent to a controlled Z rotation. This holds true 
for the other two terms. If we denote the eigenstates of 
Xby 



X\x+) 
X\x-) 



\x-) 



(41) 
(42) 



(7 ® e~ M - x ) e- 



-ia x X\ ^ia x X<S>X 

\x+}(x+\®I + \x-)(x„\®e- i2a * x , (43) 



and that 



(j (g) e ~ ia v Y ^ e iat y Y®Y 

= \y + )(y+\®I + \y„)(y_\®e- aa y Y . (44) 

These operations are equivalent to controlled rotations in 
the X and Y directions respectively. For the first case, 
if the control qubit is in the \x-) state an X rotation is 
applied to the target qubit, and not applied if the control 
qubit is in the |a;+) state. Similarly for Y. 

Single qubit rotations, Vi,V2,Wi,W2 are possible on 
the Kane quantum computing architecture, the remain- 
ing task is to specify the pulse sequence for the purely 
entangling unitary U can . Fortunately this is always pos- 
sible, as any interaction (with single qubit rotations) be- 
tween the two nuclei is sufficient [32) . In fact, it is a 
relatively simple task to use almost any interaction be- 
tween qubits to generate any desired operation. 



Calculation of the Interaction Content between 
Nuclei 



In this subsection we will see how it is possible to apply 
the canonical decomposition to the Kane quantum com- 
puter. This is important as this natural interaction of the 
system will be manipulated by single qubit unitaries to 
find the pulse scheme of any two qubit gate. The canon- 
ical decomposition provides a unique way of looking at 
this interaction. 

The interaction we will apply the canonical decomposi- 
tion to is free evolution of the configuration described in 
Sec. Ill Dl using the results cited there from second order 
perturbation theory. After a particular time of free evo- 
lution, our system will have evolved according to unitary 
dynamics, which we may decompose using the canonical 
decomposition: 



U, 



sys 



(V{ ® V 2 S ) U s can (W° ® Wi), 



(45) 



where the super-script 's' indicates a physical operation 
present in our system. 

We wish to find the interaction content U^ an of this 
free evolution. Systematic methods for doing this are 
given in |2{| |2(| |33 ■ This is most easily done by noting 
any interaction content, U can is diagonal in the so-called 
magic basis, otherwise known as the Bell basis. This 
basis is: 



6 



= ^(loo> + 


|H}), 


(46) 


— 7 

l* 2 ) = — ^(|oo> — 


|H», 


(47) 


1*3) = ^=(|oi) - 


|io», 


(48) 


|*4) = ^(|oi> + 


|io». 


(49) 


z are related to the 


eigenvalues e iAl 


e lX \ 



e zAs and e tM of U can . That is: 



Ai 

A 2 
A, 
A 4 



a 



y 



-a 

-a x + a y 
-a x - a y 
-a x + a u + a 



oe z , 
a z , 
a z , 



(50) 
(51) 
(52) 
(53) 



It is possible to relate these eigenvalues to our system. 
After a time t, each of the eigenstates of the system will 
have evolved according the Schrodingcr equation, which 
we may view as having performed an operation U sys (t) 
on the system. As we showed in Sec. Ill Dl 



where 



U svs \ll) 


= e +ieB \ll), 


(54) 


U sys \00) 


= e -^ s |oo), 


(55) 


U sys | 5) 


= e+*», 


(56) 


U sys \a) 


= e~ ies \a), 


(57) 


e s 


= u s t, 


(58) 



(59) 



Applying Eqs. - to Eqs. gHJ - 6HJ>, we obtain: 

U sys \^) = cos(0 B )|$i)-sin(0 B )|$ 2 ), (60) 

U sys \$ 2 ) = cos(0 B )|$ 2 )+sin(0 B )|$i), (61) 

U ays \$ 3 ) = e- lf,s |$ 3 ), (62) 

U sys \$4) = e+' ies |$ 4 ). (63) 

This shows that in the magic basis, U sys is given by: 



U, 



sys 



cos(e B ) sin(0 B ) 

-sin(0 B ) cos{e B ) 

e- lSs 

e 1 ^ 



(64) 



It is possible to find the eigenvalues Ai, A2, A3 and A4. 
We note that the eigenvalues of U T U in the magic basis 
are given by 



\(U T U) = {e 2lX \e 2lX \e 2lX ^e 2lXi }. 



(65) 



Calculation of the eigenvalues of Uj ys U sys is easy since 
Ug ys U S y S is already diagonal in this basis, with diagonal 
elements being {1, 1, e~ 2l9s , e 2l9s }. Care must be exer- 
cised at this point, because it is not clear which branch 
should be used when taking the argument. In our case, 
as long as < S < § [H then 



Ai = 0, 

A 2 = 0, 

A 3 = -9 S , 

A4 = +9s- 



(66) 
(67) 
(68) 
(69) 



Using Eqs. (|50|l - l|53|) we may solve for the coefficients 
a x , a y and a z , giving: 



a', 



a 



1 

-I 

2 

1 

2 
0. 



0s, 



(70) 

(71) 
(72) 



Single qubit rotations, W( , Wf , V{ , , induced are Z 
rotations. Z rotations are fast, and may be canceled in 
comparatively little time by single qubit Z rotations in 
the opposite direction: 



V^) Us 



(Wi f ® w% T ) = u t 



(73) 



For notational convenience we will now label the inter- 
action content of the system by an angle rather than by 
its time. The time for this interaction may be calculated 
through Eqs. f70)l - O, and Therefore, we 

write: 



where 



U° a n(<f>) = e^ X ® X+ ^ Y ® Y , 



(74) 



(75) 



This analysis has been based on second order pertur- 
bation theory. As we approach the electronic energy level 
crossing, this approach is no longer valid. Close to this 
crossing numerical analysis shows the eigenvalues are no 
longer symmetric which implies a s z becomes non-zero. 
Unfortunately in this regime, we excite the system into 
higher energy electronic configurations. 

Given any two qubit gate, such as the CNOT gate, 
there are many different possible choices of single qubit 
rotations and free evolution that will implement a desired 
gate. Z rotations are faster single qubit rotations than X 
and Y rotations, and therefore it is desirable to minimize 
X and Y rotations in order to optimize the time required, 
for any given two qubit gate. 
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IV. THE CNOT AND CONTROLLED Z GATES content with angles of 



A. Introduction 

The CNOT gate is a particularly often cited example 
of a two qubit gate. CNOT and single qubit rotations 
arc universal for quantum computation |34f . Many im- 
plementations, including the Kane proposal use this 
fact to demonstrate that they can, in principle, perform 
any quantum algorithm. It is a member of the so-called 
fault tolerant £3 set of gates, which are universal for 
quantum computing, and are particularly important in 
error correction. In this section we find a pulse scheme 
to implement the CNOT gate on the Kane quantum com- 
puter. 

Controlled Z rotations, sometimes known as controlled 
phase gates, are some of the most important operations 
for implementing quantum algorithms. In particular, 
one of the simplest ways to implement quantum Fourier 
transformations (QFTs ) us es multiple controlled Z rota- 
tions (see for example [36J)- Single qubit rotations and 
the controlled Z gate are, like the CNOT gate, universal 
for quantum computation. Controlled Z rotations may be 
used in the construction of controlled X and Y rotations. 
In this section we find a pulse scheme to implement any 
controlled Z rotation on the Kane quantum computer. 

Because these two gates have similar interaction con- 
tents we consider them together. Wc will first show how 
to construct a controlled Z gate of any angle, and use 
this gate directly to construct a CNOT gate. 

A controlled Z rotation of angle 9 is defined in the 
computational basis by 



U AZ (0) = 



10 
10 
10 

e ie 



(76) 



The canonical decomposition of the controlled Z rotation 
by an angle 9 has an interaction content consisting of: 



a x 
a y 

Oi z = - 



(77) 
(78) 

(79) 



This interaction content may be found by using system- 
atic methods [2£l Hi| The controlled Z gate also 
requires a Z rotation as described by Eq. (|40|l . 

CNOT is defined in the computational basis by the 
matrix 



Ucnoi 



10 
10 
1 
10 



(80) 



a x 


= o, 


(81) 


«!/ 


= o, 


(82) 




7T 


(83) 
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Since the CNOT and Controlled Z gates are both types 
of controlled rotation similar to those described in Sec. 
IIII Al it is not a surprise that they have a similar in- 
teraction content. In fact, Control Z gates (that is, a 
controlled Z rotation by an angle of tt) and CNOT gates 
have an identical interaction content, and are therefore 
equivalent up to single qubit rotations. A CNOT gate 
may be constructed from a Control Z gate conjugated by 
7® H. 



B. The Construction 

Our first task in finding a suitable pulse scheme for the 
controlled Z rotation to find a pulse scheme which imple- 
ments the interaction content [Eqs. (|77|l - (|79|) ] of the 
controlled Z rotation. Techniques have direct analogues 
in NMR [HUH!- 

The first technique is to conjugate by I®X, 7®Y ', 
or J ® Z to change the sign of two of these parameters. 
For example: 

_ g— iota-XgiX— ia. v Y®Y+ia. z Z®Z {8^) 

This can be useful, because it allows us to exactly cancel 
every controlled rotation except one: 



(I ® Z) U can (I ® Z) u c 



i2a z Z®Z 



(85) 



In our case, however, it turns out that a s z = 0. In order to 
reorder the parameters a useful technique is to conjugate 
by Hadamards |2^| . This is one of only several choices of 
single qubit rotations which reorder the parameters. In 
this case the order of the parameters is: 

(H ® H) e ia ^ x ® x + ia y Y ® Y + ia ^ z ® z (JJ g) J{~j 

_ e ia z X®X+ia y Y®Y+ia 3s Z®Z ^gg-j 

Combining these two techniques gives the following con- 
struction 

i6Z®Z 



{Z®I){H®H)U s can [-) {H®H){Z®I) 



(H ® H) U s a 



-) (H®H) 



(87) 



To find the final construction, several one qubit opti- 
mizations were made by combining adjacent single qubit 
rotations and using the identities: 



The canonical decomposition of CNOT has an interaction 



HZH 
HH 



X. 
I. 



(89) 
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Operations may be performed in parallel. For exam- 
ple, performing identical X or Y rotations on separate 
nuclei is a natural operation of the system, because mag- 
netic fields are applied globally. Performing operations in 
parallel is faster, and also higher fidelity than performing 
them one at a time. 

The construction of the controlled Z rotation is shown 
in Fig. |21 In this circuit the single qubit rotations speci- 
fied in Eq. (|4U|I have been included. The period of inter- 
action between nuclei may be increased or decreased to 
produce controlled rotations by any angle, 9, as specified 
in Eqs. (EJ), O, @ and (22). 







X 









FIG. 2: Circuit Diagram for Controlled Z Pulse Sequence 



Our task of constructing a CNOT is now comparatively 
simple. We note that a CNOT gate has the same inter- 
action term as the controlled Z (controlled phase) oper- 
ation. These gates are therefore equivalent up to local 
operations. 

Conjugation by / (£> H will turn a controlled Z oper- 
ation into a CNOT gate. Using some simple one qubit 
identities to simplify the rotations at the beginning and 
end of the pulse sequences we arrive at the decomposition 
illustrated in the circuit diagram shown in Fig. [3J 



H 








X 








H 




M^f) 




Ml) 



































FIG. 3: Circuit diagram for CNOT pulse sequence 



C. Time and Fidelity 

Throughout the paper we define fidelity as 

F(|V>),|Vo)) = |(^o}| 2 , (90) 



Description 


Time 


X rotations 
Z rotations 
2 qubit interaction 


12.6 /is 
0.2 fis 
3.2 [is 


Total 


16.0 us 



TABLE IV: Time for CNOT Gate 



with \ip) being the actual state obtained from evolution, 
and |^o) being the state which is desired. We define the 
error in terms of the fidelity as 

£ = max[l-F(|^),|V>o))]. (91) 

m 

where the maximumization is performed over the output 
of all the computational basis states, \ip). 

Numerical simulations were carried out by numerically 
integrating Schrodinger's equation for the Hamiltonian of 
the system, Eq. J2J) . The results of this numerical simu- 
lation for the pulse sequence of the CNOT gate are shown 
in Fig. 0] These graphs show each of the states and the 
transitions which are made. In these figures it is possible 
to see the evolution of each of the four computational 
basis states. The control qubit is the second qubit and 
the target qubit is the first qubit. 

According to the numerical results, a full Controlled 
Z gate takes a total time of W.lfis and has an error of 
approximately 4 x 10~ 5 . Similarly we find the CNOT 
gate takes a total time of 16.0/is. The time required for 
this gate can be grouped as shown in Table llVl 

X and Y rotations make up the majority of the time 
taken to implement the controlled Z and CNOT gates. 
In the CNOT gate, only 3.2/^s is spent implementing the 
entangling part of the gate, whereas 12.6//S is required to 
implement the X and Y rotations. 

We can see via simulation that the systematic error in 
the CNOT gate is approximately 4 x 10~ 5 . Some of this 
error will be due to errors during simulation, and break- 
down of the second order approximation. A large part 
of the error, particularly if the hyperfine interaction may 
not be varied very much, is due to X rotations where un- 
intended non-resonant transitions are excited along with 
the intended rotation. 



V. THE SWAP AND SQUARE ROOT OF SWAP 
GATES 

A. Introduction 

One of the most important gates for the Kane quan- 
tum computer is envisioned to be the swap gate. This 
is because, in the Kane proposal, only nearest neighbor 
interactions are allowed. This gate swaps the quantum 
state of two qubits. By using the swap gate it is possible 
to swap qubits until they are nearest neighbors, interact 
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2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 

Time (^s) Time (^s) 



(c) 1 10) initial state (d) |11) initial state 

FIG. 4: Numerical simulation of the CNOT gate showing different initial conditions 



them, and then swap them back again. Having an effi- 
cient method to interact qubits which are not adjacent 
to each other is therefore important, and the swap gate, 
with its high level of information transfer, is one possible 
method of achieving this. 

The square root of swap gate has been suggested for the 
quantum dot spin based quantum computer architecture 
|37| . where it is a particularly natural operation. In our 
system it is not such a natural operation, but that does 
not mean that we cannot construct it. Like the CNOT 
gate, the square root of swap (together with single qubit 
rotations) is universal for quantum computation. In this 
section we find a pulse sequence to implement both the 
swap and the square root of swap gates on the Kane 
quantum computer architecture. 



by: 



The swap gate is defined in the computational basis 



Uswap — 



10 
10 
10 
1 



(92) 



The canonical decomposition of the swap gate has an 
interaction content with angles of: 



4 

7T 

4 

7T 



(93) 
(94) 
(95) 
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FIG. 5: Circuit Diagram for the Swap gate pulse sequence 




FIG. 6: The Circuit Diagram for the Square Root of Swap Pulse Sequence 



The square root of swap gate is defined in the compu- 
tational basis by: 



Uss 



1 















4(1 +*) o 
1(1 -i) |(1 + *) o 
1 



(96) 



The canonical decomposition of the square root of swap 
gate has an interaction term consisting of: 



(97) 
(98) 
(99) 



Since the square root of swap and swap gates have es- 
sentially the same interaction content, their constructions 
are very similar, and are therefore considered together 
here. 



B. The Construction 

The easiest way to construct a swap gate is simply to 
use free evolution to obtain the angles a x and a y which 
is natural for our system. The only remaining term is 
the a z term, which for our system will naturally be 0. 
We may obtain this term by applying a pulse sequence 
similar to the Controlled Z rotation as described in Sec. 
FVl The resulting construction swap gate is shown in the 
diagram in Fig. [SJ 

The interaction content of the square root of swap gate 
is exactly half that of the swap gate, and it is negative. 



We use exactly the same technique used to obtain the 
swap gate, only allowing the nuclei to interact for exactly 
half the time. To make the terms negative we conjugate 
by Z g) I. The construction of the square root of swap 
gate obtained using this method is shown in Fig. 



C. Speed and Fidelity 

The swap and square root of swap gates were simulated 
numerically The resulting transitions for the swap gate 
are shown in Fig. {7\ Similar results were obtained for 
the square root of swap gate, not shown here. 

The swap gate takes a total time of 19.2/xs, and has a 
fidelity of approximately 7 x 10~ 5 . The majority of time 
in this gate is taken by X and Y rotations, which are also 
the major source of error. 

This is substantially faster than an existing suggestion 
for the swap gate [ijj of 192[is. It is also faster than 
using three adiabatic CNOT gates, which would take ap- 
proximately 78 (is. 

According to numerical simulation the square root of 
swap gate takes 16.8/xs and has an error of approximately 
5 x 10 -5 . This is the first explicit proposal for the Kane 
quantum computer for the square root of swap gate. 

The square root of swap gate has been suggested 
in the context of quantum computation for quantum 
dots |37| . It is universal for quantum computation 
and therefore can be used to construct a CNOT gate. 
Unfortunately in this case, a CNOT constructed from 
the square root of swap gate presented here would take 
approximately AOfis which is much longer than the pulse 
sequence presented in this paper for the CNOT gate. 
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Time (us) Time (us) 

(a) 00) initial state (b) |01) initial state 




5 10 15 5 10 15 

Time (us) Time (|is) 



(c) 1 10) initial state (d) |11) initial state 

FIG. 7: Numerical Simulation of the Swap Gate 



VI. CONCLUSION 

We have shown how the canonical decomposition may 
be applied to the Kane quantum computer. We found 
the canonical decomposition of a natural operation of the 
computer, that is. free evolution with hypcrfine interac- 
tions equal and the exchange interaction non-zero. We 
then used this interaction to form two qubit gates which 
may be applied to the Kane quantum computer. These 
gates and their times and fidelities are shown in Table Ivl 

The majority of the time required to implement each 
of these two qubit gates is used to implement single qubit 
rotations. Were we able to perform these rotations faster 
and more accurately then the gates presented here would 
also benefit. Another possible avenue of research is to 



Gate 


Time 


Error 


CNOT 


16.0/is 


4 x 10" b 


Swap 


19.2/xs 


7 x 10" 5 


Square Root of Swap 


16.2/xs 


5 x 1(T 5 


Controlled Z 


16.1/is 


4 x 10" 5 



TABLE V: Gate Times and Fidelities 



investigate the effect of decoherence on the system. 

To our knowledge this is the fastest proposal for swap, 
square root of swap, CNOT and controlled Z opera- 
tions on the Kane quantum computer architecture. We 
have shown how a representative set of two qubit gates 
may be implemented on the Kane quantum computer. 
These methods may prove particularly powerful because 
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they only involve characterization by three parameters 
which may be determined theoretically, as shown here, 
or through experiment. Once determined, these param- 
eters may be used to construct any two qubit gate. 
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